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Abstract. Exact structure function equations are an efficient means of 
obtaining asymptotic laws such as inertial range laws, as well as all measurable 
effects of inhomogeneity and anisotropy that cause deviations from such laws. 
"Exact" means that the equations are obtained from the Navier-Stokes equation 
or other hydrodynamic equations without any approximation. A pragmatic 
definition of local homogeneity lies within the exact equations because terms 
that explicitly depend on the rate of change of measurement location appear 
within the exact equations; an analogous statement is true for local stationarity. 
An exact definition of averaging operations is required for the exact equations. 
Careful derivations of several inertial range laws have appeared in the literature 
recently in the form of theorems. These theorems give the relationships of 
the energy dissipation rate to the structure function of acceleration increment 
multiplied by velocity increment and to both the trace of and the components of 
the third-order velocity structure functions. These laws are efficiently derived 
from the exact velocity structure function equations. In some respects, the 
results obtained herein differ from the previous theorems. The acceleration- 
velocity structure function is useful for obtaining the energy dissipation rate 
in particle tracking experiments provided that the effects of inhomogeneity are 
estimated by means of displacing the measurement location. 

1 Introduction 

Energy dissipation rate e is used for scaling the turbulence acceleration 
statistics that are measured by particle tracking at Cornell University 1 . The 
inertial range of the second-order velocity structure function is measured in the 
Cornell apparatus to determine e using the empirical incrtial-range relationship 
that requires local isotropy^p. In the Cornell experiments the flow between 
counter rotating blades produces high Reynolds numbers but inhomogeneous, 
anisotropic turbulence. A motivation for the present study was to use the 
exact velocity structure-function equation 2. to relate e to measurable quantities 
without making any assumptions about homogeneity or isotropy of the flow and 
hence, to obtain a relationship that applies exactly to the need to determine e in 
the Cornell experiment. By "exact" we mean that the equations follow from the 
Navier-Stokes equation and the incompressibility condition with no additional 
approximations. Because particle tracking is used at Cornell University to 
measure acceleration and velocity, it is natural to relate the acceleration- velocity 
structure function (i.e., the structure function of the product of acceleration and 
velocity increments) to both e and the third-order velocity structure function. 
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On the basis of local homogeneity, but without use of local isotropy, Mann, Ott 
and Andersen [3 E] were the first to obtain an inertial-range relationship of the 
acceleration- velocity structure function to e. By specializing our exact equation 
for the acceleration- velocity structure function to the locally homogeneous case, 
we efficiently obtain their inertial-range result in Sec. 4, but an extraneous 
derivative moment in their relationship is removed here. The exact structure 
function equation method given here not only obtains that asymptotic law, but 
also shows all quantities that must be evaluated to account for the effects of 
inhomogeneity and anisotropy. Such evaluation requires the displacement of 
the measurement volume that is defined by the Cornell laser beam and imaging 
system. 

The theorems of Nie and TanveerjS] and of Duchon and Robert [Hj es- 
tablish the inertial-range 4/3 law that relates e to the trace of the third-order 
velocity structure function. Instead of invoking local homogeneity, Nie and 
Tanveer 5 perform a large volume average, and perform a long time average 
instead of invoking local stationarity, and they use an average over orientations 
of the two measurement points instead of invoking local isotropy. Results of Nie 
and Tanveer[S] are efficiently obtained in Sec. 5 by use of the exact structure 
function equation. For the case of spatially periodic direct numerical simula- 
tion, similar results were obtained from the exact structure function equations 
in [2]. 

The theorem of Duchon and Robert differs from that of Nie and 
TanveerJSJ in an essential manner. The space-time averaging required by 
Duchon and Robert is of arbitrary extent and the viscosity is zero. In Sec. 
6, using space-time averaging and the limit of very large Reynolds number, the 
4/3 law is obtained from the exact statistical equations. The present derivation 
requires conditions on the space-time average as stated in Sec. 6 All terms that 
describe the inhomogeneity and anisotropy of flows are retained. That result 
shows all quantities that should be evaluated in direct numerical simulation 
(DNS) of flows to determine causes of deviation from the 4/3 law. 

The inertial-range 4/3 relationship of the trace of the third-order veloc- 
ity structure function is easily obtained relative to the inertial-range 4/5 and 
4/15 laws for the longitudinal and transverse components of the third-order ve- 
locity structure function, respectively. It is remarkable that the 4/5 and 4/15 
laws were obtained by Nie and Tanveerp] on the same basis as described above, 
and remarkable that the 4/5 and 4/15 laws were also proven by EyinkjJ] on a 
basis similar to that of Duchon and Robert[H], that is, for space-time averaging 
of arbitrary extent in the limit of vanishing viscosity. Those 4/5 and 4/15 
laws are not derived here. Eyink 7 finds that tests of the inertial-range power 
laws for arbitrary extent of averaging would be difficult. On the other hand, 
evaluation of all terms in an exact statistical equation by means of DNS data 
would show causes of deviations from those laws as well as the approach toward 
those laws. 

Of course, real experiments do not provide the opportunity to perform 
the above mentioned averaging used in the derivation of the theorems for the 
inertial-range laws. Deviations from 4/5 and 4/15 laws are observed because of 
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turbulence inhomogeneity, anisotropy, and finite Reynolds number. Recently, 
there has been much work |5| CH 1121 IT3*] on quantifying terms in the 
velocity structure function equations to learn how the various terms affect the 
balance of the equation as a function of scale and at what scales the effects of 
inhomogeneity and anisotropy become pronounced and how asymptotic laws are 
approached. The exact structure function equations are useful in this regard 
because they retain all effects of inhomogeneity and anisotropy in a clearly 
organized manner. The exact equations obviate the need for many derivations 
that obtain some limited aspect of inhomogeneity. Those developments are 
discussed in Sec. 7. The definition of local homogeneity that directly simplifies 
exact structure function equations is that the rate of change of statistics with 
respect to where they are measured is negligible. That pragmatic definition is 
explained and contrasted with previous definitions in Sec. 2. 

2 Pragmatic definition of local homogeneity 

Exact statistical hydrodynamics involves the derivation of equations relat- 
ing statistics without the use of approximations. The Navier-Stokes equation 
has been used to derive exact equations relating velocity structure functions 
of velocity increments and other statistics [2] • Such exact equations have been 
obtained for all orders of velocity structure functions in Structure func- 

tion equations of all orders have also been given for the asymptotic case of local 
isotropy in and ^5]. The exact equation for the scalar structure function 
has been obtained from the continuity equation as well|16j. All such exact 
equations have optimal organization when the independent variables are chosen 
to be the spatial separation r = x — x' of the two points of measurement, i.e., x 
and x', and the midpoint X = (x + x') /2. The derivative with respect to X, 
i.e., dx n i acts on some statistics within the structure function equations. Lo- 
cal homogeneity is the approximation that the rate of change of statistics with 
respect to the location of measurement may be neglected. Since that location 
is X, the result of dx n acting on a statistic is neglected relative to other terms 
in the structure function equations. This is a truly local definition; it makes no 
mention of a spatial domain. This definition of local homogeneity was exploited 
in the formulation of exact structure function equations in an d nas been 

used in studies ^3 ^2 H21 ITS*] of the effects of inhomogeneity on the balance of 
structure function equations; it was introduced in [T7|. 

Local homogeneity has been given various definitions by different au- 
thors. Kolmogorov|^j introduced a space-time domain that is small compared 
to L and T = (L/U), where L and U are "typical length and velocity for the 
flow in the whole." Kolmogorov considers the two-point differences of the veloc- 
ities at spatial points in the domain; one point is common to all the differences. 
KolmogorovJ^ defines local homogeneity as follows: the joint probability dis- 
tribution of the velocity differences is independent of the one common spatial 
point, and of the velocity at the one common point, and of time. There are 
data 1201 l2~T] that contradict the statistical independence of velocity differ- 
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ence and the velocity at either end point, and also contradict the statistical 
independence of velocity difference and the velocity at the midpoint. The 
exception is isotropic turbulence for which case local homogeneity is as- 
sured. An alternative possibility that is particularly relevant here is that the 
two-point velocity sum, u n + u' n , might be statistically independent of velocity 
difference, but statements in [21)1 I2T] contradict that statistical independence 
as well. Therefore, Kolmogorov's^^ definition should not be used because 
experimental data contradict that statistical independence 1201 EH > as do 
theoretical considerations [23 - Recently, Frisch et al. [231 have considered the 
inconsistency of Kolmogorov's definition of local homogeneity. 

Monin and Yaglom [21] define local homogeneity to mean that the joint 
probability distribution of the two-spatial-point velocity differences is unaffected 
by any translation of the spatial points. They do not impose a restriction on 
the translations to a spatial domain, ft follows that statistics composed en- 
tirely of the velocity differences obey the same relationships that they do for 
homogeneous turbulence (namely, they are independent of where they are mea- 
sured), and that the mean velocity depends linearly on positional]. In practice, 
statistics of velocity differences and of derivatives do depend on where they are 
measured except in the ideal case of homogeneous turbulence. Frisch [25] gives 
a definition that is equivalent to that of Monin and Yaglom |2I] , except that the 
translations are restricted to a domain the size of the spatial scale characteris- 
tic of the production of turbulent energy (which he calls the "integral scale"). 
Two-point structure function equations of all orders contain a statistic that is 
the product of not only factors of the velocity difference but also one factor 
of the sum of the two velocities, i.e., u n + u' n [21 E]- Because the definitions 
of local homogeneity by Monin and Yaglom ]21] and Frisch [21] involve only 
the joint probability distribution of two-point differences, but do not involve 
u n + u' n , it follows that those definitions are not sufficient to simplify structure 
function equations to the same level of simplification as does homogeneity. 

3 Definitions and notation 

3.1 Definition of two-point, two-time quantities 

Here, we extend the results in [2] to two times and to the structure function 
composed of the product of acceleration difference and velocity difference. Four 
independent variables are considered; x and x' are spatial points; t and t' are 
times. Denote velocities by = Mi(x, t), u\ = Uj(x',t'), and accelerations by 
at = flj(x, t), a[ = a,i(x',t'), etc. Quantities that will appear when we use the 
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Navier-Stokes equation are 
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where d Xi p is pressure gradient. Change of variables: 

X=(x + x')/2 and r = x - x'; T = (t + t') /2 and \ = t-t'. (1) 

Just as (x, x',t, t') are independent variables, so are (X,r,T,^). The signifi- 
cance of variable X is that it is the location of measurement. Nonzero values 
of derivatives of statistics with respect to X, evaluated at position X, are the 
result of local inhomogeneity of the flow. Likewise, nonzero values of derivatives 
of statistics with respect to T are the result of nonstationarity of the flow. The 
relationships between the spatial derivatives are 

d x% = d ri + ^d Xi , d x >. = -d Ti + 7^d Xl , d Xi = d Xi + d x >., d r% = ^ (d Xi - d x > t ) , 
which give the useful properties 

d rt [/(x, t) ± 5 (x', t')] = d Xi [/(x, t) T 9(x',t')] /2. (2) 

Similarly, 

<9 T = d t + d t >, <9 t = - (d t - d t >) . 
Use of the derivative formulas on t„ gives 

= -2(p- p') (s^ - s'ij) + d Xi [(p - p') (uj - u'j)} + d Xj [(p - p') {ui - u'i)] , 

(3) 

where = (d Xi Uj + d Xj Uj) /2 is the rate of strain. The trace of is 

Tii = 2d Xi [{p-p'){ui-u'i)]. (4) 

Use of derivative formulas on and taking the trace and use of Poisson's 
equation, i.e., d Xn d Xn p = H,n,<>, "<• gives 

e« = v- 1 [e + e') + d Xn 8 Xn (p + p'), 
where e = 2vsijSij is the energy dissipation rate. 
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3.2 Use of the Navier-Stokes equation 

The Navier-Stokes equation at (x, t) is 

di = d t Ui + u n d Xn Ui = -d Xi p + vd Xn d Xn Ui , and d Xn u n = 0. 

Subtracting the Navier-Stokes equations at (x!,t r ), multiplying by (v,j —u'j), 
and adding the equation needed to produce symmetry under interchange of i 
and j gives 

Aij = d^dij+dx n F ijn+d rn d ljn = —Tij+2v (dr n d rn dij + ^d Xn dx n d lj - e< 

(5) 

The trace gives 

An = chrdu + d Xn F «„ + d rn d lin = 2vd rn d rn du - 2 (e + e') + W, (6) 

where 

W = -2d Xi [(p - p') {ui - u'i)] + ^d Xn dx n d u - 2vd Xn dx n (p + p') ■ (7) 
No average has been used yet. 

3.3 The r-sphere and orientation averages 

Local isotropy has been used in the past to remove the divergence from 
dr n dii n and to proceed to Kolmogorov's equation. Without the assumption of 
local isotropy, an average over an r-sphere removes the divergence from d Tn dn n . 
The r-sphere average also mitigates anisotropyjJHI E]- Energy dissipation rate 
e averaged over a sphere in r-space, (&} r _ sphcl . c , was introduced by Obukhov|77| 
and Kolmogorov 28 in 1962; it is a recurrent theme in small-scale similarity 
theories. We have produced exact dynamical equations containing the sphere- 
averaged energy dissipation rate[2]. The volume average over an r-space sphere 
of radius rs of a quantity Q is defined by 

(Q)r-sphere ( X -^, T, t) = (4^1/3) ~* J J J Q (X, r,T, t) dr. (8) 

|r| < r s 

The orientation average over the surface of the r-space sphere of radius rs of 
the outward normal component of any vector q n is defined by 

£ q n = (Airriy 1 J J ^q n (X,r,t)ds (9) 

|r| = r s 

= (4^ / / ^<?„(X,r,t)dfi, 
r = rs 

where ds is the differential of surface area and dfl is the differential of solid angle. 
In JjJJl, the inner product — q n produces the so-called longitudinal component 
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of q n . The above averages can be performed on data. The divergence theorem 
is 



(<9r„gn) r _ sphcro = (3/r s ) f q n - (10) 



3.4 The X-space and X-surface averages 

Let the spatial average be over a region R in X-space. The spatial average 
of any quantity Q is denoted by (Q) H , which has argument list (R, r, T, ^); that 
is, 

(Q) R (R,r,T,t) = iy J /"QpC,r,T,li)dX, 

R 

where V is the volume of the space region R. Given any vector q n , the diver- 
gence theorem relates the volume average of dx n q n to the surface average; that 
is, 

(d Xn q n )u = 77 I I I d Xn q n dX = £ ( ~ / / iV„g„ds) = | I q n , (11) 



V J J J V \S J J J ~ v 

where S is the surface area bounding R, dS is the differential of surface area, 
and N n is the unit vector oriented outward and normal to the surface. As 
seen on the right-hand side of (|1 lf> . we adopt, for brevity, the integral-operator 
notation 

1 



x„ 



S 



N n dS. 



3.5 Time average 

Consider the term d^dij in (jSJ). The time average of a quantity Q from an 
initial time To and over a duration T is defined by 

(Q) T (X,r,T ,T,t) = - / Q(X,r,T,t) dT. (12) 

Of course, does not commute with the integral operator (|12f) ; it follows that 
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{d T d i} ) T EE - / [ftdii) dT 

•'To 

= [d y (X,r,T + T,t)-rf y (X,r,T ,t)]/r. (13) 

This shows that it is easy to evaluate (d^dij) T using experimental data because 
only the data at times T and T + T are used. One can make {d^dij) T as small 
as one desires by allowing T to be very large, provided that dij (X, r, T + T, t) 
does not differ greatly from d^ (X, r, To>t)- 
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3.6 All averages commute 

The above averages are integrations with respect to independent variables 
X, r, and T; also, an ensemble average is a sum over realizations. All those 
averages commute with one another. They commute with derivatives with 
the following exceptions: Volume and surface averages (Q) K and ^ x do not 
commute with the derivative dx t , nor do (Q) r _ sp i lcro and § r q n commute with 
d Til nor does (Q) T commute with d^. 

4 Theorem of Mann, Ott, and Andersen 

Mann et al.[3] and Ott and Mann0] studied turbulent dispersion of particles 
by means of particle tracking. They therefore used an ensemble average over 
particle trajectory events. The ensemble average is formed from observations 
of many particle-pair trajectories and r values. It is also possible to use single 
trajectories; then two times, t and t', are needed. From the data, averages are 
stored in bins of the direction and length of r. Denote the ensemble average by 
() E . The ensemble average is a sum and therefore commutes with the spatial 
and temporal derivatives. Mann et al.[5| assumed local homogeneity. From the 
present definition of local homogeneity, we have (dx n F u n ) e = {Fan) e = ®> 
and similarly (W) E — 0. Consequently, the ensemble average of © can be 
written as 

(Au) E = 8t (du) E + dr n (diin) E = 2iy dr n d rn (da) E -2(e + e') E . (14) 

For r much larger than dissipation-range scales, Mann et al. [31 neglect the term 
2vd rn d rn (du) E . Doing so, and using the definition of Ay, we have 

((a i -a' i ){u i -u' i )) E = -{e + e') E . (15) 

With two exceptions, this is the result of 3 in their equation (90) First, their 
result contains an extraneous derivative moment that is absent above. Second, 
(s + e') E appears above, whereas 2 (e) E is in their result. In an experimental 
apparatus, there can be a systematic difference between the energy dissipation 
rates at (x,i) versus (x',i'), even if t ~ t' , such that (e) E ^ (e') E . 

The above ensemble averages depend on r, which is the vector separation 
of points on trajectories from which data are obtained. That dependence on r 
allows a further average over the r-sphere. Anisotropy of the flow might favor 
the occurrence of some orientations of r; by weighting of the occurrences, the 
r-sphere average is constructed for uniform orientation of r. 

4.1 Applying the r-sphere average 

The r-sphere average causes each point x to coincide with point x' when 
the orientation of r is reversed. Consequently, 

((e + e%) h =2«e) p ) . =2((e / ) p ) , , 

x \ I hi i r-spnerc \ x 'hi' x-spnerc x \ 'hi' x -sphere 
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where averages over spheres in x-space and x'-space are introduced. The par- 
ticle tracking experiment is temporally stationary such that Dt (du) E vanishes. 
Then, the r-sphere average of 1)140 gives 

4 

\diin) E = — -Ts ((£)£) x -sphere > 

which is an extension of the 4/3 law without use of local isotropy, and 

«K " oj) («i " u'i)) E ) T _ spheie = -2 «£> B ) x _ sphere ■ (16) 

Mann et al.[3] give 1)16)1 . except that their result contains an extraneous deriva- 
tive moment. See equation (91) in 

4.2 Quantifying effects of inhomogeneity 

The terms that are neglected on the basis of local homogeneity, such as 
dx„ (F ' an) E> can be measured by particle tracking. To evaluate the effects of 
inhomogeneity, the rate of change with respect to displacement of the averaging 
volume is needed, i.e., X must be changed. The most troublesome of the terms 
describing inhomogeneity is the term —2dxi [{p — p') {ui — it^)] m because 
it requires a measurement of pressure difference. By applying the derivative 
relationship @ , the r-sphere average © , and the divergence theorem 1)1UI) , the 
most troublesome term can be expressed as 

9X< {{(P-P') {Ui -^)>B> r _ S phcro = IT f ((P+P) ( U n-U n )) E - 

'S J r n 

The right-hand side requires the two-point pressure sum correlated with the 
velocity difference and an average over orientations of r as defined inJSJ. 

5 Theorem of Nie and Tanveer 

The spatial average of (JHJ) and use of the divergence theorem 1)11)1 gives 

(^m)k = (dii) R + — <j> Fiin + d rn (da n ) m 

= 2^d rn d rn (d ll ) R -2(e + s') R + (W) R , (17) 

where 



S 



V 



x„ 



v 



(W%=77? -2(p-p')(u n -u' n )+ d Xn d ij -2vd Xn {p + p l ) . (18) 
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These averages are for any arbitrary region R. 

There are several possibilities for eliminating the terms y ^ x Fu n and 
(W) K . First, those terms describe the effects of inhomogeneity and are neg- 
ligible in the case of local homogeneity, which case requires a sufficiently large 
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Reynolds number and small |r|. Second, as in j^j, those terms vanish for the 
case of periodic boundary conditions, such as are often used in DNS, provided 
that the spatial average is over the entire spatial period. Third, those terms 
vanish if the turbulence is of limited spatial extent and the volume average is 
over a much larger region such that the velocity and pressure fields are negligible 
on the surface that bounds R; this possibility causes the statistics to decrease 
as V increases. Fourth, if the region K is enclosed by rigid boundaries with the 
no-slip boundary condition on velocities, then only the pressure gradient term 
at far right in (|18l) is nonzero at points on the boundaries; a separate hypoth- 
esis that dx n (p + p') vanishes is required. Nie and TanveerjS] eliminate 
the subject terms by integrating over the "entire volume" ; the meaning must 
be of an infinite volume unless boundary conditions are specified. Neglecting 
the subject terms, we have 

(A U ) K = <9 T {d u ) R + d rn {du„) u = 2ud rn d rn (d u ) M - 2 (e + e') R . (19) 

For r much larger than dissipation-range scales, neglect the term 2vd Tn d Tn (da) R . 
For an average over the entire flow and for t = t' , 

(e + e\ = 2{e) m = 2{e\. 

If t 7^ t' , then clearly, (e) R ^ (£')r- The definition of Aij and (|19f) give 

((o i -aO(«i-u{)> B = -2(e> R . (20) 

Apply the time average (jHj) to l(T9")l, Nie and Tanveer |S| take the 
averaging time, T, to be infinite such that (|13[1 vanishes. Alternatively, one can 
assume sufficient stationarity that the trace of (|13[1 is negligible. Also apply 
the r-sphere average (jHJ- Then (|19f) gives, for t = t' , 

f f ((^) R ) T = -^ 5 <((e} R } T ) r _ sphere . (21) 

Note that the orientation average of d iin , as defined in appears in (|2*T|) as 
a consequence of the divergence theorem (11011 . This is a theorem of Nie and 
Tanveer 5 . Their order of averaging is different, namely R then § then T, 
but the averages commute as discussed in Sec. 3.6. 



6 Theorem of Duchon and Robert 

The important distinction between the theorem of Nie and Tanveer (dis- 
cussed above) and that of Duchon and Robert[H] is that the space-time averaging 
required by Duchon and Robert is of arbitrary extent and the viscosity is zero. 
In this case, the averaging region R is a subdomain of the entire flow. Recall 
that dij = (ui — uQ (itj — u'j), and write 

Ufi — )— i/' 

F ijn — ^"nj^dij -t~ ~ dij , (^2) 



10 



where / 
U n = u n \u n )^ , u n — u n y^n/M ' ano — ^ 

Recall that (|17fl - l|18[l apply for any arbitrary space averaging region R. Let 
v = in pT fr -fTB jl . substitute (J22J, and multiply by rs/3, then 




Take the averaging region R to have a simple topology such that the volume 
to surface ratio V/S is the size of R. Because v — is effectively the limit 
of infinite Reynolds number, rg can be as small as desired. Therefore, the 
limit Srs/V — > can now be applied. This limit means that r$ is very small 
compared to the size of the averaging volume V/S, but there is no requirement 
that V/S be a length scale at which energy is generated. Thus, the size of the 
averaging volume is arbitrary provided that V / S rs- The limit Srs/V — > 
applied to f^jl gives 



+ r fd rn {d tm ) M = _|r s ^±^ . (24) 

The two terms that explicitly describe effects of inhomogeneity arc eliminated 
from H23|) . Of course, 124fl contains the following advective derivative 

Tg 1 Svg f 

Y ((un) R d x ,^) R = 3 — J x ("„) K rf« (25) 

which is seen to vanish in the limit Srs /V — > 0. The reason for not eliminating 
this advective term from l|24|l in the previous step is to make a point about 
random sweeping at the end of this section. 

The theorem of Duchon and Robert 6 applies to the case t = t'. In the 
previous section where the average was over the entire flow we had (e) R = (e') R 
if t = t' . That is not true here because R is now a subdomain of the entire flow; 
as X varies over R with r fixed, not every spatial point occupied by x coincides 
with a spatial point occupied by x' and vice versa. However, including the 
r-sphere average does cause x and x' to equally occupy every spatial point 
in the double integration. That is, for t = f ', <(£> R ) r _ sphcrc = ((e')w) r . spheIe - 
Henceforth, consider only the case t = t'. Thus, 

= (( £ )m}i-sphere = ( ( £ Wr-sphcrc ' 

r-sphcrc 
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Apply the r-sphere average (JSJ and the divergence theorem (|10|l . and neglect 
the advective term in (|25|l . and apply the time average, then (|24ll gives 

y (d T ((^) R ) r _ sphoro ) T + ^jf (diin) M ^ = -^((( £ >M) r - sp hcrc) T - ( 26 ) 



From l)13|l . the time-derivative term in l|26ll is 

y (<9 T ((rf«>R> r . sphere ) T 
rs 



3T 



«d« (X, r, T + T, 0)) R ) r _ sphoro - ((du (X, r, T , 0)} R ) r _ sphcrc (27) 



If ((<fe)]g) r _ S pi 1 ere decreases as rs decreases, then there is an averaging duration 
T sufficiently large that the left-most term in 1)26(1 may be neglected. The 
classic inertial-range case serves as an example; namely, ((du)g} r _ apileie scales 

with e 2 / 3 r^ 3 , where, for brevity, e = ((£)R) r _ sphcrc - For that case, the ratio of 
the right-hand side of (|27|l to the right-hand side of (|26[1 is not greater than 
order e _1 / 3 r^ 3 /T. That is, to neglect the time-derivative term the averaging 
duration T must be much greater than the inertial-range time scale e _1 / 3 r^ 3 . 

Recall that v = such that rs, and also e _1 / 3 r^ 3 , may be as small as desired. 
Hence T may be much smaller than integral time scales. 

Thus, with v = and a space average of arbitrary length scale, provided 
that length scale is much greater than rs, and a time average of arbitrary 
duration, provided that duration is much greater than the inertial-range time 
scale based on rs, we have 



<<Wr ) = {((e) R ) r _ sphcic ) T ■ (28) 

This is closely related to the theorem of Duchon and Robert [5]. The order of 
the averaging operations in 128|) is the same as in , although those averaging 
operations commute (see Sec. 3.6). Note that the orientation average, as 
defined in ©, appears in l|28|l as a consequence of the divergence theorem (|1U|) . 

6.1 Random sweeping 

Eyink 7 inquires into the feasibility of numerical and experimental tests 
of his theorems and the theorem of Duchon and Robert. An aspect of such 
testing is the fact that the sum of the time derivative and advective terms in 
(|24l) can be much smaller than either term taken separately. The time derivative 
and advective terms in l|24() constitute a comoving derivative. For simplicity, 
remove the factor rg/3 and the r-sphere average from these two terms; then 
their sum can be written as 

dr {du) M + {(u n } M dx n dii) R = dr {du) R + (« n ) R {dx n du) R 

= 9t {du) R + — (u„) R [ f N n d u dSj , 
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wherein two of many ways to express the sum are given. The last expression 
contains the product of the scalar da and the outward normal unit vector N n 
averaged over the surface of the averaging volume followed by the inner product 
with (u n ) R . Recall that u n = (u n + u' n ) /2; thus (it n ) R constitutes the sweeping 
velocity of the scales larger than the averaging volume as well as any mean flow 
advection. The effect of changes of the advective term is a corresponding change 
of the time derivative {du) R such that the sum of the two terms can be much 
smaller than either taken separately. 

6.2 Acceleration- velocity structure function 

The relationship of the acceleration-velocity structure function to the en- 
ergy dissipation rate is given in (|15|l and (|16|l on the basis of the ensemble 
average without and with the r-sphere average, respectively, and in (|20|l on the 
basis of an average over the entire flow. Now, (|24J) contains 

((ai-aO(ui-uJ)> R = -<E + eV (29) 

This is obtained on the basis of v = and the local space average without 
the time and r-sphere averages and without requiring t = t' . Terms describing 
random and mean- flow sweeping are absent from the derivation of 

7 Effect of large-scale inhomogeneities at small 
scales 

Recent research has quantified the influence of inhomogeneous terms on 
the balance of structure function equations. Those studies have applied to 
channel flow|l(Jl 1111 112) . wind tunnel grid-generated turbulence^ OH ^2 I12| . 
and turbulent jets^3J. The Reynolds decomposition of all terms in the exact 
equation (J5J is given in 29 . That decomposition contains all possible terms 
contributing to the effects of inhomogeneous large-scale structures. For the 
ensemble average, the Reynolds decomposition of Uj(x, t) is defined by 

u<(x,t) = Ui(x,i) +Ui(x.,t), Ui(x,t) = (ui(x,t)) E , (ui(x,t)) E = 0, 

and similarly at the point x'. Consider the Reynolds decomposition of the 
average of the term dx n Fij n in JSJ) 

dx n {Fijn) E = Un+ 2 U ' n d Xn Dij + d Xn (Aif jn + AjTj n + f iin ) , (30) 
where, for brevity, we define 

ee (dij) E , Ai = (Ui - U-) , 

Tin = ^(Ui-Uj J , 

f ijn EE ((Ui-U'i) (Uj - Uj) Un ^ ^ • 
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Note that Dij is the structure function of velocity, not of velocity fluctuations. 
For grid-generated turbulence data, one can let subscript 1 denote the down- 
stream direction and approximate U\ = U[, such that the first term on the 
right-hand side of l|3U|) is Uxdx^ij^ this is the form of the inhomogeneous term 
derived in the above-mentioned studies. |51 191 ITT1 H2] One can easily show that 
the second term on the right-hand side of i|3U|) is a generalization of the in- 
homogeneous term given in |10l ITT] for the case of uniform channel flow. Of 
course, <|3U|) and the other Reynolds decompositions given in are exact for 
every flow, whether laminar or turbulent. Equations for structure functions 
of velocity fluctuations differ significantly from equations for the velocity; that 
topic has been treated in detail in |29] . 

The exact structure functions are useful in other respects. For example, 
from 121 and (J3J it is clear that the pressure velocity term vanishes from the 
trace equations such as Ijl4|l and Ijl7|) on the basis of incompressibility and local 
homogeneity. It is not necessary to assume the more restrictive basis of local 
isotropy as in The combination of © and (|30fl shows that both dx-Dy 

and \ (U n + U' n ) dx n Dij must appear in the structure function equation as was 
correctly deduced by Lindborg [S] on the basis of mean-flow Galilean invariance, 
but replacing d^Dij with i (U n + U' n ) dx n Dij, as was done in early derivations 
on the basis of Taylor's hypothesis, does not preserve that invariance. 

8 Summary 

The formulation of the exact statistical equations in the variables X and r 
of produces the pragmatic definition of local homogeneity discussed in Sec. 
2. It is the single definition that simplifies structure function equations for 
the case of local homogeneity. The analogous definition of local stationarity 
arises from use of the temporal variables T and J of Q. The resultant or- 
ganization of the exact statistical equations provides immediate and improved 
derivation of the inertial-range relationships ljl4Jl lllfl|) between the acceleration- 
velocity structure function, the third-order velocity structure function, and the 
energy dissipation rate. Further, all terms that must be evaluated to include 
the effects of turbulence inhomogeneity and anisotropy are evident in the for- 
mulation. That makes the exact statistical equation method useful to particle 
tracking experiments such as the turbulence acceleration experiments at Cornell 
University. The recent theorems of Nie and Tanveer |S] in l|21ll and of Duchon 
and Robert [B] in l|28l) are likewise obtained easily, as are their generalizations, 
e.g., H17[). The relationship between the acceleration- velocity structure func- 
tion and energy dissipation rate is obtained in (|29() on the basis of only a local 
spatial average for v = 0. With (|30J) as one example of the Reynolds decom- 
position of all terms in the exact statistical equations in |29| . it is evident that 
the exact statistical equation formulation efficiently reveals all terms that must 
be evaluated to quantify the influence of inhomogeneity on the balance of struc- 
ture function equations. It is not necessary to derive the individual terms that 
describe the effects of inhomogeneity that are missing from equations valid only 



14 



for homogeneous turbulence. All such terms are now known. 

The above are examples of the use of exact statistical equations. Their 
usefulness arises from the organization of the equations. 
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